The role of pressure on the magnetism of bilayer graphene 
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We study the effect of pressure on the locahzed magnetic moments induced by vacancies in bilayer 
graphene in the presence of topological defects breaking the bipartite nature of the lattice. By using 
a mean-field Hubbard model we address the two inequivalent types of vacancies that appear in the 
Bernal stacking bilayer graphene. We find that by applying pressure in the direction perpendicular to 
the layers the critical value of the Hubbard interaction needed to polarize the system decreases. The 
effect is particularly enhanced for one type of vacancies, and admits straightforward generalization 
to multilayer graphene in Bernal stacking and graphite. The present results clearly demonstrate that 
the magnetic behavior of multilayer graphene can be affected by mechanical transverse deformation. 

PACS numbers: 



I. INTRODUCTION 

The magnetic properties of graphene and its multilayer 
compounds remain one of the most interesting topics in 
the system that still awaits experimental confirmation. 
The improved experimental capabilities to produce and 
manipulate large samples has given rise to a renewed in- 
terest on the issuei It is known that vacancies, edges 
and other defects that lead to dangling bonds in the 
graphene system induce localized magnetic moments that 
might give rise to interesting applications. The theoret- 
ical paradigm around the magnetism in graphene is the 
Lieb theorem^ that fixes the spin of the ground state of 
the bipartite system to be half the number of the un- 
paired atoms in the lattice. Although the theorem is 
demonstrated only for the Hubbard model in bipartite 
lattices, the result has proven to be very robust and to 
hold in more general calculations based on ab initio or 
density functional methods^i^. While the interactions 
can be extended, it has been recently proven that the 
ground state magnetization is very sensitive to the pres- 
ence of local topological defects such as five or seven 
rings breaking the bipartite character of the lattice^^i^ 
In particular it has been shown that, when one of the va- 
cancies inducing magnetic moments is reconstructed to 
form a pentagon, the critical value of the Hubbard in- 
teraction needed to reach a finite polarization increases 
significantly.^ 

Bilayer graphene (BLG) is even more interesting than 
single layer graphene (SLG) under many points of view,^ 
in particular for the magnetic properties. In the Bernal- 
stacking, BLG can support two types of vacancies giving 
rise to unpaired atoms: these produced by removing a 
site having a neighbor in the adjacent layer are named 
/3, and these coming from sites that are not connected to 
the other layer called a vacancies. In a bipartite lattice, 
unpaired atoms give rise to zero energy states and the 
physics underlying the magnetic properties of the system 
is that of the electronic interactions in the manifold of 
zero energy states. It has been recently shown^'^^ that 
the significant differences between the wave functions of 
the zero modes associated to the a and /3 vacancies in 



BLG give rise to different physical behaviors. In partic- 
ular when the system is gapped by an external gate, the 
vacancies of type a generate fully localized states inside 
the gap. 

Recent experimental progress in production and ma- 
nipulation of graphene samples have broadened the possi- 
bilities of tailoring the properties of SLG and BLG. Pres- 
sure is known to play an important role on the stability 
of the gate induced gap in BLG^^ and on the impurity 
stated. 

In this work we analyze the effect of pressure on the 
behavior of the localized magnetic moments coming from 
the two types of vacancies in the presence of a topolog- 
ical defect. We show that the critical interaction value 
Uc needed to polarize the system decreases for increasing 
pressure. The effect is particularly enhanced for vacan- 
cies of type a. We clarify the physical mechanism leading 
to this behavior and propose pressure as a way to improve 
the magnetism of the samples. 

We work with the tight-binding (TB) model, including 
a Hubbard interaction. The results are obtained by a self- 
consistent computation within the unrestricted Hartree- 
Fock approximation. A further analytical analysis using 
degenerate first order perturbation theory explains the 
results in terms of the different energy lifting that occurs 
in the subspace of the zero modes under the perturba- 
tions originated by the pentagonal hopping. We argue 
that the physical result will remain when going beyond 
the simple calculation performed in this work and extend 
the results to other graphene multilayers. 



II. MODEL 

The lattice structure of AB Bernal-stacking BLG is 
schematically represented in Fig. [IJa). The atoms of 
the A sublattice in the top layer are connected by t± 
to these of the B sublattice of the bottom layer. We 
will consider the minimal model for the BLG with only 
intralayer t and inter layer t± couplings. The estimated 
values of these parameters in the system are t ~ 3eV, 
and ti_ ^ 0.3 - 0.4 eV - 
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Figure 1: (Color online) (a) Bilayer lattice structure and main 
tight-binding parameters. Left: a vacancy. Right: {5 vacancy, 
(b) Reconstructed vacancies modeled by a pentagonal link. 



The two different types of vacancies are also shown in 
Fig. [1] (a). The a vacancy where the removed atom in 
the top layer {Bi in the figure) is not connected to the 
second layer is represented in the left hand side. The 
right hand side represents a {3 vacancy {Ai in the fig- 
ure). Fig. [1] (b) shows the reconstructed pentagonal 
link named tp in this work. 

The non-interacting TB Hamiltonian for the tt- 
electrons of the bilayer system isJ^ — 
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HtB = Y.^^-^^Y. [4a(R)^2a(R) + h.c], (1) 
i=l R,cr 

where Hi corresponds to the SLG Hamiltonian: 

Hi = 

-t al{-R) [6,,(R)+6,,(R-ai)+6,,(R-a2)] +h.c. , 

R,cr 

(2) 

and aicr(R) [6icr(R)] are the annihilation operators for 
electrons at position R in the sublattice Ai (Bi) of the 
layer i (i = 1,2) with spin polarization cr. The basis 
vectors can be chosen as ai = aCx and a2 = (i{e^ — 
A/3ey)/2, with a = 0.246 nm being the lattice spacing. 

The lattice of BLG is a bipartite lattice and the TB 
minimal model has electron-hole symmetry. Vacancies 
are modeled by suppressing the corresponding lattice site. 
No reconstruction or relaxation will be included in the re- 
maining structure which maintains its original geometry. 
This approximation does not affect the zero-modes we 
are interested in^^. We will consider the simplest topo- 
logical defect modeled by first producing a vacancy and 
then adding a pentagonal link connecting two of the clos- 
est atoms to the vacancy as shown in Fig. [TJb). Since 
the two atoms belong to the same sublattice this defect 
breaks the bipartite nature of the lattice. 

The magnetic behavior of BLG in the presence of va- 
cancies, edges, and other defects has been investigated 
usine; the Hubbard model in The interacting 



TB Hamiltonian is H = Ht b + Hu , with 

Hu = U Y[na,^{R)nac^{R) + nb,t{R)nui{R)] , (3) 

R,L 

where nxia{R) = ^L(^)^^cr(^), with x = a,6, ^ = 1,2 
and a =t, |. 

We use finite clusters with periodic boundary condi- 
tions at half- filling (one electron per atom). The Hamil- 
tonian is solved in the Hartree-Fock approximation, and 
the mean-field spin density at each lattice site is ob- 
tained self-consistently. It is known that a finite stag- 
gered magnetization appears in the honeycomb lattice 
above a critical value of the on-site Coulomb interaction 
Uc ~ 2.2&^. We will keep the values of U below this 
value to make sure that the physics explored is due to 
the magnetic moments associated to the defects on the 
lattice. Evidence of the presence of localized magnetic 
moments around isolated vacancies has been reported re- 
cently in scanning tunneling microscopy experiments on 
a graphite surface^^. 



III. VACANCY-INDUCED ZERO-ENERGY 
STATES 

In the honeycomb lattice, a vacancy gives rise to a 
quasilocalized state with a continuum limit wave function 



where K and are the reciprocal space vectors of the 
two inequivalent corners of the first Brillouin zone, and 
(x, y) are distances from the vacancy position^^ 

An analytic expression for the vacancy-induced states 
in the continuum model of BLG was obtained recently 
in Ref. 10 following the procedure outlined in|To.24 for 
the monolayer case. By cutting the lattice into left and 
right regions with respect to the vacancy position a zigzag 
edge to the left and a Klein edge to the right appear. 
The wave function is obtained by matching surface state 
solutions at the zigzag edge with that localized at the 
Klein edge. Depending on the type of vacancy, two so- 
lutions are obtained. The f3 vacancy produces a zero- 
energy state quasilocalized on atoms of the opposite sub- 
lattice in the same layer of the vacancy, decaying as 1/r 
away from the vacancy. Assuming the vacancy to be lo- 
cated in layer 1, the continuum limit wave function can 
be written as 

Ti{x,y) ^ "^{x.y) 

r2{x,y) ^0, (5) 

where specifies the wave function component on 
layer and "^{x^y) is the quasilocalized state given in 
Eq. The zero-energy states induced by a vacancy of 
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OL type has the wave function 
Ti(x,^) ^{x,y) 

T2(x,7/) ^ ^e-^2V^•'^ + ^e^2^e^^'•^ (6) 

where d = arctan(?//x). This is a delocahzed state, with 
the pecuharity of being quasilocahzed in one layer (where 
the vacancy sits) and delocahzed in the other We no- 
tice here that the wave function of the /3 vacancy is in- 
sensitive to while this parameter enters explicitly in 
the wave function of the a vacancy, a fact that will be 
important in the forthcoming analysis. 

IV. NUMERICAL RESULTS 
A. Magnetism and the pentagonal links 

Since the lattice of the BLG model used in this work 
is bipartite the magnetic behaviour of the system follows 
the Lieb's theorem as in SLG. Therefore in the presence 
of vacancies, for a repulsive value of the Hubbard interac- 
tion U the ground state of the system at half filling has a 
total spin equal to half the number of unbalanced atoms 
5* = (N A — Nb)/"^' We have analysed the magnetic be- 
haviour of the two different types of vacancies that occur 
in the BLG in the presence of a pentagonal link. 

We observed that both types of vacancies have the 
same behaviour and it coincides with what happens in 
SLG^. We have considered two vacancies on the same 
layer and reconstructed one of them by forming a pen- 
tagonal link tp as pictorially shown in Fig. [IJb). In this 
situation a finite critical value of the on-site Coulomb in- 
teraction is needed to reach the ground state polarization 
predicted by Lieb's theorem. For values of U below the 
Uc the total spin of the ground state is zero. By varying 
the pentagonal hopping integral from O.lt to t we have 
verified that the critical value of U increases monotoni- 
cally with tp for both types of vacancies. 

B. The role of pressure 

We have explored the influence of pressure on the mag- 
netic properties of BLG. The effect of increasing pressure 
is modelled by increasing the hopping between atoms of 
different layers t^. The most common value for con- 
sidered in BLG is ~ O.lt although recent works con- 
sidering different types of stacking in few layers graphene 
have produced values up to t±_ ^ 0.502 eV with t = 3.16 
eV^. 

In order to exemplify the behaviour found in this work 
and for simplicity we will study the magnetic polariza- 
tion of the system in the situation described before: Two 
vacancies of the same type, one of them reconstructed 
forming a pentagonal ring. We have seen that the effect 
of pressure is very different for the two different types of 
vacancies of BLG. 
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Figure 2: (Color online) Critical value of the Hubbard pa- 
rameter U separating zero from finite polarization regions as 
a function of the strength of the pentagonal linktp, for differ- 
ent values of the perpendicular hopping t± in units of t. The 
dots are computed values; lines are a guide to the eye. 



We have varied the perpendicular hopping in the range 
t± = O.lt to t± = t and for each value of t± we increased 
the pentagonal hopping tp from O.lt to t. We found that 
the magnetic behavior of the system is independent on 
t±. The value of the Uc increases with tp but does not 
depend on t±. Pressure has no effects on the spin polar- 
ization of the system in this case. This can be understood 
by observing that, since the wave functions induced by 
/3 vacancies have amplitude only on one layer the behav- 
ior is just that of the SLG. This behavior will change if 
next nearest interlayer hoppings 73 and 74 are taken into 
account, with /3 vacancies becoming more similar to a 
vacancies. 

Considering a vacancies we observed that for a fixed 
value of the pentagonal hopping tp, Uc decreases as t± 
increases. This indicates that pressure could help the po- 
larization of the BLG ground state. The critical values of 
the on-site Coulomb interaction are represented against 
the pentagonal hopping in Fig. [2l Each curve corre- 
sponds to a different value of t±. While the strength of 
the pentagonal hopping conspires against the polariza- 
tion, the perpendicular hopping favors the appearance of 
magnetic moments. In the next section Sec.[Vlwe explain 
this behavior doing an analysis with degenerate pertur- 
bation theory. 

We have also considered case of two vacancies in differ- 
ent layers. To obtain a non zero magnetization according 
to the Lieb's theorem the two vacancies must belong to 
the same sublattice, A or B. Therefore in this case to get 
a magnetic ground state both types of vacancies a and 
P must be present. We observe that t± has no effect on 
the critical value of the Hubbard interaction needed to 
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obtain the spin polarization of the ground state when the 
pentagonal link is attached to a vacancy of j3 type. On 
the contrary, when the pentagonal ring is formed by re- 
constructing the a vacancy Uc decreases as increases 
as it occurs in the case of having both vacancies of a 
type. Therefore the increase of affects the magnetic 
behavior when a pentagonal hopping exists. Otherwise 
the Lieb's theorem holds for any value of t^. 



V. PERTURBATIVE ANALYSIS 

In the following we explain the results obtained in the 
previous section by studying how the zero-energy modes 
introduced in Sect [TTIl are affected by the pentagonal link 
giving rise to the topological defect shown in Fig. [T]^b). 
We use first order degenerate perturbation theory. The 
discussion is put on general grounds having in mind ap- 
plication to multilayer graphene and graphite. 



Definitions 



B. First order degenerate perturbation theory 

We want to know what happens to zero energy modes 
once we add the perturbation defined by Eq. (p!Q|) to the 
system. Here we use a perturbative analysis, and apply 
degenerate first order perturbation theory within the zero 
energy mode sector. This is justified, strictly speaking, 
if the perturbing parameter tp is much smaller than the 
HOMO-LUMO gap. 

The new energies for |?/;* ) within first order perturba- 
tion theory are given by the eigenvectors of the matrix 



T 



/Tn Ti2 ... 



^22 



(11) 



V 



where 



T,, = {^i\H,\^i) = 







7 G A 



(12) 



Let H be the non-interacting TB model Hamiltonian 
for some graphitic system preserving the bipartite nature 
of the lattice, as is the case of Eq ([T]). Let as assume the 
system to hold n vacancies belonging to the same sub- 
lattice, with n induced zero energy modes |?/;^) . . . 
such that 



if |v;) = o. 



(7) 



Explicitly, using localized atomic orbitals as TB basis, 
{|1) , . . . , |A^)}, where N is the number of lattice sites in 
the system, we can write the zero energy modes as 



(8) 



as we have assumed the vacancies to belong to the A sub- 
lattice. 

Owing to the quasilocalized nature of vacancy induced 
zero modes (in multilayer systems we assume to act 
on the layer where the zero mode has its quasilocalized 
component) we may consider the limiting case where va- 
cancies are sufficiently far apart. In this particular case 
zero modes will be almost unaffected if tp connects two 
sites that are also sufficiently far apart from any vacancy. 
Once these two sites approach a given vacancy we will see 
the energy of the associated zero mode go up as 



(13) 



Without loss of generality, we assume the vacancies to 
belong to the A sublattice, which allows us to write 



(9) 



jeB 



In the case of SLG, for example, the continuum limit 
would give l?/^* ) ^ ^(x,^), with ^{x^y) as given in 
Eq. dH). 

We want to perturb the system with a term that dis- 
rupts the bipartite nature of the lattice. The simplest 
such term is just an extra local hopping connecting the 
same sublattice. 



(10) 



bital I7), and ^ is a vector connecting next nearest neigh- 
bors. The local perturbation defined by Eq. (p!Q|) provides 
a simple parametrization of the pentagonal link shown in 

Fig.HTb). 



where and are the real amplitudes of that par- 

ticular zero mode at the perturbed sites. The other zero 
modes being almost unaffected. This is certainly a good 
description for the pentagonal link shown in Fig. mb), 
when tp connects two of the closest sites to a given va- 
cancy. 

It was shown in Ref. that in SLG the pentagonal link 
induces a finite Hubbard interaction Uc to polarize the 
system. Within the present perturbative analysis this 
can be understood as a consequence of the energy shift 
of the zero mode affected by the pentagonal link: the 
effect of the Hubbard term has to overcome the energy 
scale set by the shift in Eq. (fT3|) . 

Multilayer graphene, and in particular BLG, are of es- 
pecial interest. There we can tune by applying pres- 
sure, since by increasing hopping between layers the am- 
plitude of the zero mode over the layer where the vacancy 
resides decreases. This is apparent in Eq. (|6|) for zero 
mode induced by the a type vacancy and explains the 
numerical results shown in Fig. [2l 
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VI. SUMMARY AND DISCUSSION 

The possibility of tuning magnetic behavior by lattice 
deformation in graphitic materials is nowadays being con- 
sidered. Tunable magnetism by mechanical control in 
graphene is a hot topic due to the broad field of ap- 
plications of organic magnets and, from the conceptual 
point of view a challenge. In this work we have shown 
that pressure applied in the perpendicular direction to 
the planes can override the negative effect of topological 
defects breaking the sublattice symmetry in multilayer 
graphene. 

Vacancies are the principal type of defects produced in 
graphitic materials by ion bombardmen t^^i^^ and play a 
substantial role in material properties. Topological de- 
fects breaking the sublattice symmetry are also energet- 
ically favorable and are being considered in defect engi- 
neered devices^S — . 



We have seen that the two different types of vacan- 
cies that can form in Bernal stacked multilayer graphene 
have different magnetic behaviors under pressure. Both 
of them are affected by topological defects in the sense 
that a larger value of the Hubbard interaction is needed 
to polarize the ground state. We have seen that the criti- 
cal U decreases by applying pressure to the sample in the 
case of having a majority of a vacancies. The /3 vacancies 
are not affected by pressure. These different behaviors 
are due to the differences in the wave functions of the 
zero-energy states induced by vacancies of the two types. 
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